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Z '• Abstract 

, As a follow-up on our earlier paper we calculate the x^ of various AW po- 

p^ I tential models with respect to the np scattering data. We find that only the 

^P ■ most recent potential models give a reasonably good description of these data. 

Almost none of the potentials is truly an AW potential in the sense that it 
gives a good description of both pp and np scattering simultaneously. 
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In a recent paper |l|, we investigated the quality of a number of AW potential models 
with respect to the pp scattering data below 350 MeV. By comparing to the pp data we 
tested the isovector {1 = 1) partial waves of these models. In this Brief Report we will 
confront the models with the np data so that we can also test the isoscalar (/ = 0) partial 
waves. 

In our previous comparison [0, we restricted ourselves to a comparison with respect to 
the pp data only, because of two reasons. Firstly, at the time of our comparison we were 
still in the process of analyzing the np data, so we did not have a complete and consistent 
np database. Secondly, our initial goal for comparing NN models was to bring attention to 
the fact that these models are often not AW models at all, but rather pp or np models; i.e., 
they were only fitted to one type {pp or np) of data. Especially in the last decade, it has 
been assumed by many people that it is sufficient to fit the parameters of a potential model 
only to the np data. Simply adding the proper electromagnetic interaction is then assumed 
to provide the corresponding pp version of that potential, so that one in fact ends up with 
an AW potential. (This is much easier than fitting to the pp data directly, because in fitting 
only the np data one does not have to worry about all kinds of electromagnetic corrections 
when calculating the observables.) For example, the full Bonn potential was originally 
presented as an NN potential, although its parameters were in fact fitted only to the np 
data. In a later publication Q], the model had to be adjusted to make it also applicable 
to pp scattering. We wanted to show explicitly that in order to construct an WW potential 
model it is generally not good enough to only fit the np data. 

One often argues incorrectly that the reason for this failure is only artificial, where the 
argument goes as follows. There are some very accurate pp differential cross sections |Q and 
np total cross sections [^ at low energies, which put very tight constraints on the pp and np 
^Sq scattering lengths. It is well known that the nuclear pp and np scattering lengths are 
different and that this difference cannot be explained by mass differences alone. It is then 
obvious that the choice of fitting a potential to either the pp or np scattering length will 
result in a totally incorrect np or pp scattering length, respectively, for that model. Very 
high x^ values on the pp data for models which are fitted to the np scattering length will 
then be the result. 

A way to avoid this problem is to exclude these low-energy data and only include data 
above 5 MeV, say. Fitting to either the pp or np scattering length causes only minor 
differences in the quality of the description of the data at energies above 5 MeV. However, 
in doing this we found |]T| that most of the potential models still gave a very poor description 
of the pp data. This indicates that the poor quality of these models is not only due to an 
incorrect scattering length, but has other sources as well. 

Another argument that a high x^ on the pp data is only artificial is that one can slightly 
adjust one of the parameters of the potential model to ensure that it fits the pp ^Sq scattering 
length. Because the scattering length is extremely sensitive to small changes in the potential, 
one expects that such an adjustment will be minimal and will hardly affect the other partial 
waves. This modification is then believed to improve the description of the low-energy data 
without affecting the description of the other data, resulting in an overall good description 
of the pp data. We investigated this [|I| extensively using the Arg84 potential 0. Instead 
of modifying one of the parameters, we completely replaced the pp ^So phase shift of the 
Arg84 potential by the pp ^Sq phase shift of the Nijmegen partial-wave analysis 0. This 



corresponds to having an almost "perfect" ^5*0 phase shift; something which is very unhkely 
to be achievable by simply adjusting one of the potential parameters. Although this gives 
a considerable improvement, the resulting x^/^data is still rather high. The same trend is 
found when we apply this procedure to the coordinate-space version of the Bonn potential 
or to the coordinate-space Bonn A and Bonn B potentials ^j. This indicates that the poor 
quality of these potential models with respect to the pp data is not only due to the fact 
that they have an incorrect ^5*0 phase shift. In fact, because of the mass difference between 
neutral and charged pions the pp and np isovector phase shifts in the other partial waves 
are significantly different as well (see Table I of Ref. ^). Therefore, the problem usually 
cannot be fixed by simply adjusting one of the parameters of the model. 

Now that the Nijmegen partial- wave analysis of the np data below 350 MeV is finished [§ , 
we have a complete and consistent np database, and we can complete our comparison of 
NN potential models by confronting them with the np data. This allows us to investigate 
whether np potential models which give a very poor description of the pp data can in fact 
survive a quality test with respect to the np data. As we will see, some of them do. This once 
more strongly supports our statement that "a good fit to the np data does not automatically 
guarantee a good fit to the pp data" [|l|. 

In this Brief Report we will confront a number of potential models with the np scattering 
data. This means that we have to calculate the phase shifts of the lower partial waves at all 
the energies at which experimental data are available, construct the scattering amplitude, 
and calculate the observables. Since we are only going up to Tiab = 350 MeV, it is sufficient 
to calculate the phase shifts due to the nuclear potential up to total angular momentum 
J = 6. For the higher partial waves we can take the phase shifts as given by one-pion 
exchange. The np database we use is given in Ref. 0, with the difference that here we will 
leave out the low-energy data below Tiab = 5 MeV as discussed earlier. The database now 
consists of 2458 data between 5 and 350 MeV. 

The potential models we will consider are mainly the same as those of our pp compar- 
ison 0, namely: HJ62 fig, Reid68 |l|], TRS75 pl, ParisSO 0, UrbSl H, Arg84 p, 
BonnR [§], Bonn89 @J^, and Nijm93 [|l^. The soft-core Nijmegen potential is here referred 
to as Nijm93. This is an updated version of the Nijm78 potential |I^ . 

We made a change in notation with respect to the Bonn potentials, in that the coordinate- 
space Bonn potential is now denoted by BonnR, and Bonn89 refers to the full Bonn potential. 
Because we do not have the computer codes to calculate the phase shifts of the full Bonn 
potential ourselves, these phase shifts were obtained from the computer software SAID [|T7 . 



We take the phase shifts in steps of 1 MeV as generated by the SAID program, and then 
use linear interpolation to obtain the phase shifts at the experimental energies. Incidentally, 
this provides another reason for only using the np data above 5 MeV: The phase shifts below 
10 MeV in the SAID program are represented by an effective-range parametrization which, 
especially for the S waves, is not good enough to represent the low-energy phase shifts to a 
high accuracy. We checked for the other potentials that starting at 5 MeV rather than at 
10 MeV does not make any difference for the conclusions we will draw below. 

As a check that the way how we treat the full Bonn potential (Bonn89) is justified, we 
applied the same procedure to the Nijmegen potential: We took the phase shifts of the 
Nijm93 potential in steps of 1 MeV as generated by the SAID program, interpolated to get 
the phase shifts at the experimental energies, calculated the total x^ with respect to the 



np data, and compared to the result we obtain when we calculate the Nijni93 phase shifts 
directly using our own computer code for the Nijm93 potential. The difference in x^ is less 
than 0.3. This strongly supports our belief that for the present purpose we can use the phase 
shifts of the Bonn89 potential as stored in the SAID program. 

The results for the nine potential models are presented in Table |. Similarly to our 
comparison of potential models with respect to the pp scattering data [|ll , only a few of the 
models give a reasonable description of the np scattering data. The best ones are the Arg84 
and Nijm93 models with x^/-^data ~ 2, followed by the UrbSl and Bonn89 models with a 
still reasonable x^/^data ~ 3. To get some insight in where the high y^ for the other models 
comes from, we divide the total x^ in a set of sub-x^ for each different type of observable, 
as listed in Table |. 

The relatively high x^('^tot) for most models is mainly due to one group of total cross 
sections |jl8|, consisting of 70 data from 39 to 350 MeV. The statistical error on these data is 



rather small and the energy range is very large. So if the energy dependence (the shape) for 
the total cross sections as given by a potential model is different from what it is implied to be 
according to the experiment, the x^ rapidly increases. Indeed, the other total cross-section 
experiments cover much smaller energy ranges and the description of these data is generally 
much better with x^/^data ranging from 1.1 for the Arg84 potential to 3.7 for the BonnR 
potential. However, there still remains the fact that this group of 70 data can be described 
very well in the Nijmegen partial- wave analysis with x^ = 56.8. 

High x^(o"(6')) and x^{^yyi^zz) indicate an incorrect ^Pi phase shift and ei mixing pa- 
rameter. For the HJ62, Reid68, and Urb81 models the ^Pi phase shift becomes too negative 
at high energies, whereas for the BonnR model it does not become negative enough. At 
^lab = 300 MeV the Nijmegen partial-wave analysis gives 5{}Pi) = —27.58(22)°. Sim- 
ilarly, comparing to |^ ei = 4.03(17)°, the ei mixing parameter is too high in the HJ62, 
Reid68, and TRS75 models, whereas it is too small in the Bonn models. 

As an example of the differences between the phase shifts of the various potential models, 
we show in Fig. |l]the ei mixing parameter and the ^D2 phase shift. The solid curves are the 
values from the multi-energy Nijmegen partial- wave analysis 0. The shaded band represents 
the statistical error. Compared to the Nijmegen analysis, the Paris80 potential gives values 
which are too high, whereas the Bonn89 potential gives values which are too low. This 
is reflected in the relatively high x^ on the differential cross sections and spin- correlation 
paremeters for these models. 

In Table | we also give the x^/^data for the 5-350 MeV range and 0-350 MeV energy 
ranges. The latter is included to clearly demonstrate the enormous changes that can occur 
when we include the low-energy data. Obviously, including these low-energy data in a 
potential comparison can be misleading in that the large total x^ obscures the fact that a 
certain potential model can in fact describe the scattering data at energies above 5 MeV, 
say, reasonably well. 

To summarize, we find that only the Arg84 and Nijm93 potentials give a good description 
of the np data (x^/Adata ~ 2), while the Urb81 and Bonn89 potentials can be qualified 
as reasonable (x^/Adata ~ 3). When we include the very accurate low-energy total cross 
sections, only the Nijm93 potential still has a good x^/-^data- 

Before we can make statements about the general quality of A7V potential models, we 
have to recapitulate the quality of these models with respect to the pp scattering data. 



Rather than referring to our previous pubhcation [||], we here present the results for the pp 
data using the same procedure we use for the np data. So also for the pp data, we here give 
X^/iVdata results from a direct comparison to the pp data, rather than using the Nijmegen 
representation of the x^ hypersurface of the pp data. The results are given in Table |I[ Note 
once more the enormous rise for some of the models (Arg84 and BonnR) when we include 
the low-energy (0-5 MeV) data. The quality with respect to the pp data for each model 
individually has already been discussed in our previous publication [Q]. 

Comparing Tables H and |I|, we conclude that most potential models are not NN models 
at all, but give only a reasonable description of either the pp or the np scattering data (or 
not even that). The only exception is the Nijm93 potential. The situation for the Bonn89 
potential (except at very low energies) is probably also still not unsatisfactory, although the 
description of the np data is considerably worse (but not too bad) than the description of 
the pp data. The TRS75 potential also shows a quality which is consistent for both pp and 
np scattering, although x^/-^data = 3.5 on all data is not too good. All other potentials 
can be classified as either pp potentials (Reid68, ParisSO), np potentials (UrbSl, Arg84), or 
they do not give a satisfactory description of either pp or np scattering (HJ62, BonnR). Two 
other coordinate-space versions of the Bonn potentials, Bonn A and Bonn B [||, also belong 
to the last category. For the 5-350 MeV energy range, x^/^data for these two models is: for 
Bonn A, 9.4 on pp and 8.3 on np, and for Bonn B, 8.5 on pp and 8.9 on np. It is surprising 
to see the enormous difference in quality of the full Bonn potential (Bonn89) on the one 
hand and all the Bonn coordinate-space versions on the other hand. The difference cannot 
be explained p|,H by claiming that coordinate-space potentials are necessarily of inferior 
quality. The quality of the Nijmegen potential, which is a coordinate-space potential (but 
which also has an exactly equivalent momentum-space version), clearly contradicts this. 

The results from Tables | and ^ are rather disappointing, considering the fact that all 
models were originally presented as being AW potentials. This demonstrates once again 
our point that in general one has to be careful when using these potential models in other 
calculations, like in few-nucleon scattering and bound-state calculations, pp bremsstrahlung, 
or nuclear matter calculations. In most cases these models cannot even describe both pp 
and np scattering with the same, satisfactory, quality. 

We should mention that recently there have been constructed a number of new NN 
potentials which are truly NN potentials in the sense that they give an excellent description 
of both the pp and np data simultaneously. These are two Nijmegen potentials Nijm I and 
Nijm II [^], a regularized update of the old Reid68 potential |TB|, and an update of the old 
Arg84 potential [|19[. All four models have the almost optimal x^/-^data ~ 1 on both pp and 



np data. Part of the success of these models is that they explicitly contain the one-pion- 
exchange potential with the proper neutral- and charged-pion masses. Another reason is of 
course that these models were explicitly fitted to both pp and np data simultaneously giving 
the proper constraints on the isovector partial waves. 

The work of V. Stoks was supported by the Australian Research Council. V. Stoks is 
also grateful to I.R. Afnan for discussions and comments on the manuscript. 
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TABLES 

TABLE L x^ on the np scattering data between 5 and 350 MeV for various AW potential 
models. A division is made showing the sub-^^ on the total cross sections, the differential cross 
sections, the analyzing powers, the spin-correlation parameters, the depolarizations, and the rota- 
tion parameters. The lower part shows the x^/^data for the 5-350 MeV and 0-350 MeV energy 
ranges. 





A'data 


HJ62 


Reid68 


TRS75 


ParisSO 


Urb81 


Arg84 


BonnR 


Bonn89 


Nijm93 


O"tot 


225 


2376 


6599 


389 


1117 


594 


662 


1624 


895 


554 


aie) 


1323 


3379 


5771 


3273 


2856 


2726 


2055 


14908 


3138 


1554 


Ay 


738 


2345 


4217 


3799 


1275 


1807 


2019 


3330 


1791 


1023 


yy ' -^zz 


86 


975 


9456 


1054 


3899 


1331 


231 


5822 


1328 


1499 


Dt 


43 


56 


133 


87 


108 


81 


51 


296 


90 


52 


At,Rt 


43 


60 


127 


176 


202 


145 


82 


348 


89 


100 


all data 


2458 


9190 


26303 


8778 


9457 


6684 


5100 


26328 


7330 


4783 


5-350 MeV 


2458 


3.7 


10.7 


3.6 


3.8 


2.7 


2.1 


10.7 


3.0 


1.9 


0-350 MeV 


2514 


4100 


3740 


2925 


3300 


3400 


3.3 


10.5 


423 


1.9 



TABLE II. x^/^data on the pp scattering data for the 5-350 MeV and 0-350 MeV energy 
ranges. 



5-350 MeV 
0-350 MeV 



^data 

1590 

1787 



HJ62 Reid68 TRS75 ParisSO Urb81 Arg84 BonnR Bonn89 Nijm93 



9.7 
13.5 



2.5 
2.9 



3.3 

3.4 



2.2 

4.5 



5.9 
6.0 



6.9 
7615 



12.4 
1090 



1.8 
25.5 



1.9 

1.8 



FIGURES 

FIG. 1. ei mixing parameter and ^D2 phase shift of the Nijmegen partial-wave analysis with 
the statistical error (shaded band), the Paris potential (dotted line), the full Bonn potential 
(dash-dotted line), and the Nijmegen potential (dashed line). 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9411002vl 
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